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While galaxies move through the intracluster medium of their host cluster, they experience a ram pressure which removes 
at least a significant part of their interstellar medium. This ram pressure stripping appears to be especially important for 
spiral galaxies: this scenario is a good candidate to explain the differences observed between cluster spirals in the nearby 
universe and their field counterparts. Thus, ram pressure stripping of disk galaxies in clusters has been studied intensively 
during the last decade. I review advances made in this area, concentrating on theoretical work, but continuously comparing 
to observations. 
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1 Introduction 

Galaxies populate different environments, ranging from the 
field, characterised by a low galaxy number volume density, 
to densely populated clusters. In nearby clusters, especially 
cluster spiral galaxies differ from their counterparts in the 
field in a number of properties (see also the extensive review 
by Boselli & Gavazzi|2006 and references therein): 

- cluster galaxies are HI deficient compared to their field 
counterparts. The deficiency increases towards the clus- 
ter centre. Spatially resolved studies reveal that the HI 
deficiency is caused by a truncation of the gas disks. 
While the HI disks of field spirals typically extend be- 
yond the optical disks, the opposite is true for HI defi- 
cient spirals. 

- Luminous cluster spirals have, on average, a lower star 
formation rate (SFR). The suppression of star formation 
(SF) goes hand in hand with HI deficiency. Spatially re- 
solved studies reveal that e.g. the Ha disks are also trun- 
cated. 

- Cluster spirals tend to be redder than field spirals, also 
indicating less active SF. 

- An increased 20 cm radio continuum intensity suggests 
an increase in magnetic field (MF) strength by factor of 
2-3 compared to field spirals. 

- late-type galaxies follow more radial orbits and tend to 
have higher velocites than early-type galaxies, which 
suggests they are free-falling into cluster 

The global contents of molecular gas seems to be compara- 
ble between field and cluster spirals, the investigation of the 
global dust contents is difficult and still controversal. 

All taken together, these observations suggests that one 
or more processes in cluster environments remove gas from 



galaxies or make them consume their gas, which leads to 
a subsequent decrease of SF activity and hence change in 
colour. A number of processes are suspected to be responsi- 
1980| l suggested the starvation scenario: 



ble: Larson et al. 



With a moderate SFR, a typical spiral galaxy consumes its 
disk gas within a few Gyr. Dense environments cut the disk 
galaxies off its external gas supply, e.g. their halos, and thus 
SF ceases. Tidal interactions or mergers are more common 
in galaxy group environments, and have a profound impact 
on the affected galaxies. In clusters, encounters between 
galaxies happen at a much higher velocity, thus the interac- 
tion time is much shorter. Nonetheless, [Moore et al.| ( |1996[ 
1 1998) 1 showed that harrassment, the repeated high velocity 
encounters combined with the tidal field of the cluster, can 
affect cluster galaxies in the desired way. Finally, galaxy 
groups and clusters do not only contain galaxies, but also the 
intragroup or intracluster medium (ICM). Especially cluster 
galaxies pass through this medium at high velocity, typi- 
cally slightly supersonic ( Faltenbacher & Diemand|[2006 



and thus experience a substantial ram pressure (RP) on their 
gas disks (or gas halos). By comparing the average gravita- 
tional restoring force working on the gas disk to the RP ex- 
pected in massive clusters like Coma, Gunn & Gott (1972, 
hereafter GG72) estimated that typical spirals should have 
their gas pushed out by the ICM head wind. Thus, the idea 
of ram pressure stripping (RPS) was bom. 

Naturally, RPS is relevant not only for disk galaxies in 
clusters, but in a variety of contexts, e.g. for elliptical galax- 
ies (e.g. [Lucero et al.[2005t[McCarthy et al.|2008| l and dwarf 
galaxies in clusters and groups, and even dwarf galaxies in 
the gaseous halos of giant galaxies (e.g. the Large Magel- 



lanic Cloud, M astropietro et al.[[20()8t |Mastropietro|[2009 
pLghi et al.|2009 
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By comparing the different processes, [Boselli & Ga-[ 
|vazzi[ ([2006 ) concluded that RPS is the best candidate to 
explain the observed differences between field and cluster 
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spirals in the nearby universe. As shown below, RPS is able 
to remove gas from a galaxy in the desired way, the SFR 
and colour change in the aftermath. Besides being able to 
explain the trends in global properties, several galaxies are 
known to show detailed characteristics expected for RPS 
candidates (Sect.|4|. 

RPS of disk galaxies in clusters was subject to intensive 
studies during the last decade. In this review, I will focus 
on recent work dedicated to RPS of massive disk galaxies 
in galaxy clusters. The emphasis is put on theoretical work, 
though I will frequently compare the theoretical results to 
observations. I will start by summarising the results of basic 
theoretical work in Sec. |2] Subsequent sections deal with 
different extentions to the basics in the effort to arrive at a 
coherent picture. Sec.[TO]will summarise the current status. 

2 Basic theory: the gas removal process 

2.1 The Gunn&Gott criterion for tlie stripping radius 

In the simplest picture, the amount of gas lost from a gala- 
xy's disk is determined by the competing forces of RP and 
gravitational restoring force. The case where the disk galaxy 
moves face-on through the ICM is easiest to access: Con- 
sider a disk galaxy in cylindrical coordinates, {R, Z). In or- 
der to determine the restoring force for a given disk radius, 
R, one needs to find the gravitational acceleration in Z- 
direction, az, i.e. perpendicular to the disk. For a symmetri- 
cal potential, az directly in the disk plane is zero. Thus, for 
a given R, one can move up from the disk in Z-direction and 
find the maximum of az for this R. Multiplying this max- 
imum gravitational restoring acceleration, azma.x{R), with 
the gas disk's surface density, yields the graviational 

restoring force per unit area. 



{R) ■ 



(1) 



at this disk radius R, which can now be compared to the 
RP. The restoring force per unit area is a decreasing func- 
tion of disk radius. For a given RP, the radius at which the 
restoring force equals the RP is called stripping radius. In- 
side this radius, the galaxy's gravitation outweighs the RP 
and thus can prevent gas removal, outside this radius the 
RP wins and strips the gas away. This method for estimat- 
ing the stripping radius is commonly called Gunn&Gott cri- 
terion (G&G criterion, for short), although it contains two 
important improvements compared to the original version 
of GG72: First, it is spatially resolved instead of bein aver- 
aged over the whole galaxy. Second, the original calculation 
of the gravitational restoring force used the mass of the stel- 
lar disk only, whereas the method described above includes 
gravitational forces due to all components of the galaxy, 
i.e. disk, bulge and dark matter halo. Like the original ver- 
sion, it is restricted to galaxies moving face-on through the 
ICM. 

The RP 



is determined by the density of the cluster gas, picM, and 
the galaxy's velocity through this medium, Wgai. Thus, RPS 
is expected to be most effective near cluster centres, where 
both, picM and Wgai, are highest. 

2.2 Early simulations 

The first RPS simulations were dedicated to the removal 



of gas replenished by stars from elliptical galaxies (IGisler 
[TW; 'Lea & Young" 1976'; Shaviv & Salpeter"1982' Takeda 
et al. 1984; Gaetz et al. 1987, Portnoy et al..l993, Balsara 
et al.|19941 i 



The first simulations of disk galaxies (Farouki & Sha 
|piro||1980{ [Toyama & Ikeuchi||l980, ) investigated whether 
RPS can convert gas rich spirals into gas poor SO. Despite 
low resolution and the restriction to two dimensions, already 
these simulations agreed that for typical cluster conditions 
RPS can strip typical spirals severely or completely, and that 
stripping works at the outer parts of the disk first. 

2.3 Systematic studies 



Later work (|Abadi et al.|1999t Schulz & Struck 


2001 Mar- 


colini et al. |2003[ Roediger & Hensler 2005 


Roediger & 


Briiggen 2006") performed more systematic studies at higher 



resolution, both, in 2D and 3D. Simulations were done mostly 



for massive disk galaxies, but also for dwarfs (MarcoUni 
|et al.|[2003] ), and mostly for ICM conditions representative 
for cluster centres, but also for low RP environments ( |Roedi 



ger & Hensler|2005j Roediger & Briiggen 2006 , Marcohni 



et al.|2003| l. The key question was to test the G&G criterion 



for the stripping radius. Generally, for face-on geometries, 
this simple criterion does a good job. 

Studying galaxies inclined w.r.t. the orbit requires 3D 
simulations. Test cases were included by Abadi et al. ( 1999|l; 
IQuilis et al.|(|2000|;|Schurz & Struck|(|2001|l.|Marcolini et al ' 



j2003^ performed a systematic study for dwarfs, 'Roediger 
|& Briiggen (2006 )for giants. The common result is that in- 
chnation does not have a significant influence unless the 
galaxy moves close to edge-on. Mass loss in face-on strip- 
ping (defined as inclination 0°) up to an inclination of 60° is 
very similar. The inclination is relevant mostly for medium 
RPs, that strip a galaxy severely, but leave some gas in the 
inner part of the disk. Strong RPs that strip a face-on galaxy 
completely, also strip an edge-on moving galaxy, although 
on a somewhat longer timescale. Weak RPs that hardly af- 
fect a face-on galaxy, neither strip an edge-on galaxy. The 
most prominent influence of inclination is that stripping of 



inclined galaxies makes the gas disk asymmetric (Roediger 
l& Bruggen.2006) . 

2.4 Timescales: ram pressure stripping as a 
multistage process 



More recent simulations fSchulz & Struck'2001', 'Marcolini 



etal. 



PlCM^^gal! 



(2) 



2003[ |Roediger~& Hensler 2005 ; Roediger & Briiggen 



2006 1 which were run for longer than 100 Myr, revealed that 
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RPS is actually a multi-stage process. The first stage could 
be called RP pushing: here the RP pushes the disk gas from 
disk radii larger than the stripping radius. This displacement 
happens on short timescales of a few 10 Myr 

However, the gravitational potential of a galaxy is deep 
and extends beyond the disk region. Consequently, acceler- 
ating the displaced disk gas to escape velocity requires an 
ongoing RP. Unbinding the disk gas from the galactic po- 
tential takes a few 100 Myr 

In addition to the gas removal due to the RP, grid sim- 
ulations can resolve the Kelvin-Helmholtz (KH) instability 
that is expected to work continuously on the gas disk's sur- 
face. Gas loss due to this continuous or turbulent-viscous 
stripping ( |Nulsen||1982 1 happens at a much lower rate (< 
IMq yr~^) and thus becomes evident in long-time simula- 
tions only. 

2.5 Common conventions and simplifications 

The "classic" RPS as proposed by GG72 happens only dur- 
ing the first of the phases described in the previous section, 
it is the same as the RP pushing. In recent years, also the 
complete multistage process is referred to as RPS. Alterna- 
tively, the complete process is called gas stripping or ICM- 
ISM interaction (ISM for interstellar medium). 

Simulations are mostly done in the galaxy rest frame, 
where the motion of the galaxy through the cluster translates 
into an ICM wind flowing past the galaxy. Therefore, the 
equation for the RP often reads picMVicw where the ICM 
wind velocity, wicm. replaces the galaxy's orbital velocity. 

Different methods are used to model the ICM-ISM inter- 
action. The simplest are particle codes, where the RP is in- 
cluded as an additional force on gas particles. Hydrodynam- 
ical simulations utilise either smoothed particle hydrody- 
namics (SPH) or grid codes. Each method has its strengths 
and problems (e.g. Agertz et al.|2007 Tasker et al.|2008[ l. 

Although the RP varies along the orbit through a clus- 
ter, the first step was to use constant ICM wind. This sim- 
plification was used by all work discussed in this Sec. |2] 
The second important simplification is that of a homoge- 
neous gas disk, although the interstellar medium (ISM) of 
real disk galaxies is highly inhomogeneous. On the obser- 
vational side, the first indication regarding the gas contents 
was the observation that many cluster spirals have truncated 
gas disks. The homogeneous gas disks in the simulations are 
commonly meant to be comparable to the diffuse HI disks. 

The next sections describe different extentions to this 
basic work. 

3 Varying ram pressure during cluster 
passage 

When a galaxy passes through a galaxy cluster, both, its or- 
bital speed and the ICM density surrounding it, increase to- 
wards the cluster centre and decrease after centre passage. 
Consequently, the temporal RP profile is peaked. 



An obvious step to improve the basic simulations is to 
make the ICM wind vary accordingly. For elliptical galax- 
ies, this was attem pted since the 70ies (|Lea & Young|I976[ 



Takeda et al.|1984[[Ac"reman et al.|2003| l, although with re- 
strictions to 2D, low resolution and radial orbits. |Toniazzo| 
& Schindler ( 200I[ l modelled the passage of an elUptical 



galaxy on a rosetta-like orbit through a Coma-like cluster in 
2.5D. 

The constant ICM wind simulations for disk galaxies 
suggest interesting effects. Given that it takes a few 100 
Myr to accelerate the gas displaced in the RP pushing phase 
enough to become unbound from the galaxy's potential, a 
short enough RP peak may be unable to unbind all gas, but 
a fraction of the displaced gas could fall back to the disk. 
This would mean a deviation from the G&G criterion. 

Such a fallback indeed takes place in the simulations of 
[Vollmer et al.| ( |200T] l, who were the first to model RPS of a 
disk galaxy in a varying RP. In this work, the gas disk was 
modelled by a sticky-particle code. The RP was included 
as an additional acceleration on particles at the upstream 
side of the galaxy. This description cannnot capture hydro- 
dynamic effects like underpressure in the wake, eddies, or 
KH instabilities. The orbits through the cluster were approx- 
imated to be purely radial, i.e. the RP does not vary in direc- 
tion, allowing to use of an analytical temporal RP profile. 
The model cluster is tailored to resemble the Virgo clus- 
ter. As Virgo is very compact, the RP peak is rather short. 
In this case, stripping becomes a distinct event, and indeed 
displaced gas falls back to the disk when the RP ceases. The 
authors propose that this feature is useable to distinguish be- 
tween ongoing and recent stripping. 

Roediger & Briiggen (2007, hereafter RB07) and [JF 
|chym et aL] ( |2007| ) were the first to study the cluster passage 
with full 3D hydrodynamic simulations. The SPH simula- 
tions of Jachym et al. ( 2007[ l; Jachym et al. ( 2009) concen- 
trated on galaxies orbiting on purely radial orbits through 
rather compact, Virgo-like clusters and thus also found RPS 
to be a distinct event, leading to backfall of displaced gas. 

The grid simulations of RB07 were dedicated to a com- 
pact cluster (comparable to Virgo, though not as compact) 
and an extended. Coma-like cluster. These simulations were 
conducted in the cluster rest frame. While most of the clus- 
ter is covered in a low resolution grid, the adaptive grid al- 
ways refines highly on the moving galaxy (see Fig. [TJ. In 
the galaxy's wake, the resolution was kept on the 1-2 kpc 
level, allowing a study of RPS tails (see Sec.|5]l. Only at dis- 
tances to the galaxy larger than 100 kpc, derefinement was 
enforced to limit computational costs. Despite the lower res- 
olution at larger distances from the galaxy, this method even 
aUows to track the distribution of the stripped gas through 
the cluster In this more general case studied in these simu- 
lations, no backfall of gas specifically after core passage is 
seen. Only in case of a very rapid increase of RP, the gas 
loss is somewhat delayed compared to the expectation of a 
time-dependent G&G criterion. One reason may be that the 
cluster studied by Vollmer et al. is more compact than the 
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Fig. 1 RPS simulation of RB07,RB08. Slice through the 
computational grid in the orbital plane, showing the colour- 
coded gas density. The background picture shows the whole 
cluster, the inset zooms on the galaxy. White boxes show 
blocks of 8'^ grid cells. The galaxy was resolved with 250 
pc. 



ones studied by RB07, another may well be the differing 
numerical methods. 

A comparison of the simulated mass loss history with a 
time-dependent version of the G&G criterion revealed that 
this simple criterion is a reasonable approximation in the 
obvious situation: while the galaxy falls into the cluster, the 
galaxy should not move near edge-on, and the RP should not 
increase too rapidly. Otherwise simulations and analytical 
estimate differ (Fig. [2]). Moreover, the mass loss due to KH 
stripping in long low RP or edge-on periods is not captured 
by the G&G criterion. 

Considering RP histories for various orbits in the two 
clusters reveals interesting implications for the cluster ga- 
laxy population (RB07): The constant ICM wind simula- 
tions distinguished strong (> 10" 
10^^^ ergcm^'^), and weak (< 10" 
gimes, where a typical massive spiral (rotation velocity 200 
km s~^) would be stripped completely, seriously (say down 
to a few kpc gas disk radius) or just marginally. In realistic 
clusters, no orbits with only weak RPs are possible. Even 
in the compact cluster on orbits with large impact param- 
eters, a galaxy experiences at least a medium RP. In order 
to move on a more circular orbit, a galaxy needs to have a 
higher orbital veloctity, which results in a medium RP de- 
spite a low ICM density. For radial orbits, the high central 
ICM density together with the high central velocity nearly 
always leads to complete stripping. Especially in the ex- 
tended. Coma-like cluster, typical spirals are stripped com- 
pletely before even reaching the cluster centre. Briiggen & 



erg cm medium (~ 
"i^ergcm^'^) RP re- 
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Fig. 2 Evolution of the gas disk radius (top, thick lines) 
and the gas mass (middle) during a cluster passage. The bot- 
tom panel displays the evolution of the true inclination, that 
is, the angle between the galaxys rotation axis and the di- 
rection of motion. The red and blue lines are for two galax- 
ies with different inclination. The thick black line is the es- 
timate from the G&G criterion. For the galaxy moving at 
medium to face-on inclination, the G&G criterion describes 
the mass loss history well, for the edge-on galaxy it does 
not. The thin lines in the top panel show the maximum and 
minimum radius of the gas disk at each time. The edge-on 
galaxy is asymmetric, the other is not. From RB07. 

4 Studies on individual galaxies 

4.1 Observable signatures of ram pressure stripping 

The main impact of RPS on a disk galaxy is the trunca- 
tion of its gas disk down to the stripping radius, whereas 
the stellar and dark matter components are left completely 
unaffected. In contrast, gravitational interactions tend to in- 
fluence all components of the involved galaxies. Thus, an 
undisturbed stellar disk combined with a truncated gas disk 
is commonly interpreted as a sign that this galaxy has suf- 
fered RPS. Additionally, the stripped gas should form a tail 
on the downstream side of the galaxy and thus make it pos- 
sible to determine the galaxy's direction of motion. Rotation 
curves and 2D velocity fields of the remaining gas disks are 
expected to be only weakly effected by RPS, and only in 
near edge-on cases ( [Kronberger et al.|2008b[ ). 



De Lucia ( 2008 1 investigated RP histories of cluster galax- 4 2 Examples 



ies in more detail by analysing the Millenium Simulation 
( Springel et al.|2005| ), finding that about one quarter of ga- 
laxies in massive clusters are currently subject to strong RP, 
and more than 64% of current cluster members have expe- 
rienced strong RPs since they fell into their host cluster. 



Especially Vollmer et al. (1999 - now, for a summary see 
|VoIImer|2009j ) have done extensive work in observing can- 
didate galaxies deeply in HI, radio and other wavelengths. 
Additionally, they modelled each candiate with a sticky par- 
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tide code (Vollmer et al. 2001[ l, aiming at disentangling 
each galaxy's PRS history and stage. The authors stress that 
this aim cannot be reached by comparing projected gas den- 
sity maps alone, but velocity information is essential. Al- 
though this method is applied successfully to a number of 
galaxies, one needs to keep in mind that the applied numer- 
ical method lacks some essential physics. E.g. the fallback 
of not fully stripped gas after core passage is used as a cri- 
terion to distinguish ongoing and recent RPS, the use of this 
criterion is not supported by hydrodynamical simulations 
(RB07). 

The (incomplete) list of examples given here contains 
members of Virgo, the most nearby and hence best-studied 
cluster This list focusses on the remaining (gas) disks, RPS 
tails are discussed in Sec. |5] 



NGC4522 is one of the best-studied RPS candidates, data 
in many wavebands is available. All wavebands detecting 
ISM ( HI: [Kenney et al.|[2004] l. Ha: [Kenney & Koopmann 
[T9991 IVollmer et al.||2000[ CO: |W lmer et al.| |2008| (po- 
larised) radio continuum jVollmer et al. 2004bTfar infrared: 



upstream edge (Abramson & Kenney'2009'). From simula- 
tions, ^\^llmerj (^2009^) suggests pre-peak stripping. No sig- 
nificant enhancement of SF is found, and only very little SF 
in the tail. 



NGC4569 was studied extensively by Boselli et aL (|2006 1 
and others (e.g. |Wilson et al.|2009||Chyzy et al.|2006[|Vollm4 r 
et al. 2004a| . Also for this galaxy data at many wavelengths 



is available. From multi-wavelength data, |Boselh etal. (2006 1 
date the removal of the outer gas disk to to have happened 
100 to a few 100 Myr, in agreement with the dynamical 
model of lVoUmer et aL] ( |2004a| l. 

5 Ram pressure tails 

Whereas the physical picture for the gas removal appears 
to have converged, the fate of the stripped gas is less clear, 
both from theoretical and observational point of view. 

5.1 Theoretical work 



Murphy et aL]|2009[ |Wong & Ken ney 2009 ) consistently First theoretical attempts were again made for elliptical ga- 



show a gas disk truncated at 3 kpc disk radius, whereas the 
stellar disk is mainly undisturbed (Fig.|3]l. Moreover, signif- 
icant fractions of extraplanar Ha and HI are found, suggest- 
ing a tail of 3kpc length. The polarised radio continuum sug- 
gests compression of the ISM at the upstream edge. From 
optical and UV data, Crowl & Kenney ( 2006[ l show that the 
SFjust outside the stripping radius ceased just ~ 100 Myr 
ago, after a modest star burst. 

Comparing sticky particleH-MF advection (see Sec. |7]i 
simulations to these observations, |VolImer et al.| ( |2006] l con- 
clude that this galaxy currently suffers RPS, or has experi- 
enced the RP maximum just 50 Myr ago. 

Thus, NGC4522 appears to be a textbook example of 
RPS. However, the needed RP to produce the observed strip- 
ping is unusually high given the large cluster-centric dis- 
tance (800 kpc) of this galaxy. Kenney et al. ( 2004[ i and 
Vollmer et al. ( 2006) 1 suggest that local enhancements in 



laxies ( [Stevens et al.|1999|[Tbniazzo & Schindler|2001||Acr^ 



ICM density or bulk flows may be responsible for the en- 
hanced RP at this galaxy's position, a plausible explana- 
tion, because the Virgo cluster is still in the stage of assem- 



|man et aL 2003| l, predicting weak, but observable X-ray 
wakes. The basic studies on disk galaxies naturally see the 
stripped gas forming a tail on the galaxy's downstream side. 

Studies aiming at the morphology of the tails, however, 
need to consider a number of issues generally not included 
in the work discussed so far: (i) Studying the tail requires 
a larger volume in full 3D to be simulated with sufficient 
resolution, which causes higher computational expenses for 
both, SPH and grid simulations, and thus was not done in 
most cases, (ii) The tail structure will depend on the RP his- 
tory, (iii) The tail structure will depend on additional physi- 
cal properties of ICM and ISM, like viscosity, thermal con- 
duction, cooling, self-gravity, SF, and MFs. Approaches so 
far have realised only subsets of these requirements. 

The first study dedicated to RP wakes ( jRoediger et al.| 
|2006| l found flaring tails of stripped gas, which must how- 
ever be interpreted very carefully, as a constant ICM wind 



was used. The subsequent simulations (RB07, [Roediger^fe 



bly. Investigations of cosmological simulations (Tonnesen 



& Bryan|2008t [Ludlow et al.|200"9| ) support this idea by re 



Briiggen|2008b hereafter RB08) modelled the passage of a 



vealing that at a given clustercentric distance, galaxies can 
experience RPs of a range up to or more than one order of 
magnitude, and that cluster galaxies can have unorthodox 
orbits. 

NGC4402 shows a truncated HI and radio continuum 
disk. Moreover, ablation of dense clouds seems to be at 
work (see |Crowl et al.|2005| and references therein). 



NGC4430 shows truncated HI ( [Chung et al.p007| . Ha, 
near UV, IR and radio continuum, and radio deficit at the 



galaxy though a cluster in a 3D adaptive grid code (Fig.[T]l. 
Assuming that the simulated gas disk is HI gas and remains 
so after stripping, mock HI maps can be produced (Fig.|4]i. 
For a sensitivity limit of ~ 10^^ cm"^ in projected gas 
density, typical tail lengths are 40 kpc. Such long tails are 
seen even at large distances (0.5 to 1 Mpc) from the cluster 
centre. Tail length and density are not necessarily largest 
in the cluster centre, where the RPs is strongest. Higher 
RP causes a higher gas loss per time. However, higher RPs 
are also accompanied by high velocities. Consequently, the 
stripped gas is distributed over a larger volume. The mass 
loss per orbital length, which depends on RP and orbital 
velocity, determines density and (observable) length of the 
tail. 



www.an-joumal.org 



© 2009 WILEY- VCH Verlag GmbH & Co. KGaA, Weinhdm 



6 



E. Roediger: Ram pressure stripping in galaxy clusters 



mc 4522 
HlonR 







CO (solid) 




J. v HI (dotted) 

CI 







RIGHT ASCENSION [J2000) 



RIGHT ASCENSION g20()01 



Fig. 3 Textbook example NGC4522: an undisturbed stellar disk, but truncated HI disk and extraplanar gas. The Ha and 
CO disk are truncated consistently with HI. Stellar disk: R-band, greyscale in left panel. HI: contours in left and middle, 



dotted right panel. Ha: greyscale middle panel. CO: solid contours in right panel. (Left and middle from Kenney et al. 
20041 right from |Vollmer et al.|2008| l. 
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Fig. 4 Simulated RPS tails (RB08): projected gas density 
contours. Contour spacing is logarithmic, half an order of 
magnitude per step. The cyan coloured (the 4th from out- 
side) contour corresponds to ~ 2 • 10^^ cm^. The black line 
marks the orbit. 



at distances between 0.6 and 1 Mpc from the cluster centre. 
Comparable cases were found in the simulations of RB08, 
which showed a similar structure in projected gas density, 
but generally the simulated wakes show a much more tur- 
bulent wake than observed. At first glance, long tails at such 
large cluster-centric distances are surprising, because RPS 
is strongest near the clusters centre. However, as discussed 
above, the observability of the tail depends on the gas loss 
per orbital length, not RP alone. As also discussed above, 
galaxies at a given cluster-centric radius can experience a 
whole range of RPs, so these examples could be fast-moving 
galaxies. 



Using an SPH code including gas cooling and recipes 



for SF, but a constant ICM wind, |Kapferer et al.| ( |2009| l also 
produce mock observations. Their simulated tails tend to be 
very long, up to a few 100 kpc. Part of the gas in the tail 
continues to cool and form dense knots and stars all along 
the tail. 

Although the current simulations seem to be able to re- 
produce aspects of the observations (see below), they should 
be interpreted cautiously. E.g. the physics of SF itself is not 
yet fully understood, and the relevance of MFs and transport 
processes remains to be investigated. 



Near N GC4388, a Virgo spiral, [Oosterloo &" van Gor-| 
kom 2005 observe a 120 kpc long HI cloud that appears 



to be, at first glance, a RPS tail. However, neither the tail 
structure in the HI image nor the velocity structure along 
the tail can be explained by the current simulations (RB08). 
Whether this means that this tail is not RP induced or the 
simulations are insufficient, remains to be seen. 



5.2 Comparison to Observations 

Also the observational situation remains unclear Given the 
many galaxies with truncated HI disks, at least some HI 
in tails or simply distributed through the clusters could be 
expected. However, this gas is not found easily. |VolImer| 
& Huchtmeier (2007) searched the vicinity of a number 
of known RPS candidates deeply and found no more HI 
than already known. A blind HI survey (j Kent et aL]|2007 



In addition to HI tails, there are only very few galax- 



Gavazzi et al.|2008] l found only a handful HI clouds unas- 
sociated with galaxies. One of them could be interpreted as 
RP stripped gas ( Kent et al. 2009| l. A deep targeted survey of 
Virgo spirals ( Chung et al.|2007 1 found HI tails for galaxies 



ies that have X-ray (Wang et al. 


2004 |Sun et al.|2006|l. Ha 


(Gavazzi et al. 2001 


Sun et a 


.|2007a| Yagi et al. 2007 1 


or radio tails (e.g. Gavazzi et al. 1995). These tails tend 



to be rather long (several 10 kpc), and extremely narrow 
(< 10 kpc) and straight. Moreover, especially the X-ray tails 
seem to be rare: Sun et al. ( |2007b| l have searched for addi- 
tional cases in the Chandra and XMM data of 62 galaxy 
clusters and did not find any. The few known X-ray and 
Ha tails are generally much narrower and much straighter 
than the simulated tails. Thus, additional physics like a vis- 
cous ICM, the influence of cooling and tidal effects may be 
needed to explain the observations. 
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6 Ram pressure stripping of a multiphase 
interstellar medium and the effect on star 
formation 



Answering this question observationally is rather dif- 
ficuh, as for no galaxy the true velocity w.r.t. the ICM is 
known. 



It is commonly argued that RP stripped galaxies will also 
cease to form stars and thus subsequently become dead and 
red, because the galaxy has lost its gas reservoir, which 
would be needed for further SF activity. However, the sim- 
ulations discussed so far mostly dealt with the diffuse gas 
phases of the ISM (diffuse HI, ionised) only. A realistic ISM 
is highly inhomogeneous and contains also HI and molecu- 
lar clouds. As stars form in molecular clouds, the link from 
loss of diffuse gas to SF quenching has to be made very 
carefully via a multiphase ISM. 

Modelling RPS of a multiphase gas disk poses quite a 
challenge. A variety of physical processes needs to be in- 
cluded: heating and cooling of the ISM, self-gravity, SF and 
stellar feedback, MFs. An enormous range of spatial scales 
needs to be resolved, starting from subparsecs to resolve 
molecular clouds, the sites for SF, to at least 100 kpc, the 
size of RP tails. Given that the process of SF itself is still 
a matter of ongoing research (e.g. Bate|2009 Price & Bate 
2009 Krumholz et al.|[2009 1, and given limits for compu- 
tational expenses, current simulations can only be approxi- 
mations. 

The question of RPS of a multiphase ISM has two sides: 
how does RPS differ between a homogeneous gas disk and 
an inhomogeneous one, and how is the multiphase ISM and 
thus the galaxy's SFR affected by RPS. 



6.2 Effect of ram pressure stripping on star formation 

All observations and simulations agree that at outer disk 
radii, where the diffuse gas is stripped, also SF ceases. Al- 
though the stars form in dense clouds, the internal dynamics 
of the ISM connect the cloudy and diffuse phase and thus SF 
to the diffuse phase. 

For the remaining disk, the answer is not so easy. The 
first-glance expectation is that the high ICM pressure leads 
to compression of the ISM and thus increased SF. Earlier 
work ( |Fujita & Nagas hima 1999] |Bekki & Couch|2003l ) in- 
dicates an increase in SFR by a factor of a few. A similar 



or even stronger effect is seen by Kronberger et al. ( 2008a i 
and |Kapferer et al.| ( [2008||2009| , who utihse an SPH code 
including cooling and standard recipes for SF and stellar 
feedback. Interestingly, stars are also formed in the stripped 
gas. However, their model galaxies experience a constant 
RP. Thus, the evolution of the SFR in a varying RP is still 
an open question. 

Observations do not give a clear answer, either. A study 
on Ha properties of Virgo spirals ( Koopmann & Kenney] 
2004c ) Koopmann et al. 2006 Koopmann & Kenney 2006 ) 



suggests that the SF in the outer disks, where gas was stripped, 
ceases, whereas the SFR in the remaining inner disk is mod- 
erately enhanced. However, other observations do not find 
an increased SF despite clear RPS signatures (NGC4330 
Abramson & Kenney|2009| l. 



6.1 Gas loss from a multiphase disk 

The first-glance expectation for the first question is that an 
inhomogeneous disk will be stripped more easily, because 
the ICM wind blows away the low surface density parts of 
the disk, creating holes. This in turn enlarges the surface of 
the disk, which is more vulnerable to KH instabilities. RPS 
should be more efficient. However, in the light of non-linear 
internal dynamics of ISM, one could argue towards an oppo- 
site answer: The increasing ICM pressure (staticH-RP) may 
compress ISM clouds, increasing surface density, making 
them harder to strip. Also MFs may prevent KH instabili- 
ties. A first attempt to find an answer by simulations was 
done by [Quilis et"ar ( 2000 1, who artificially cut holes into 
an otherwise smooth gas disk. No coohng was included, 
thus the disk is stripped more easily. A more advanced at- 
tempt was made by |Tonne sen & Bryan] ( |2009| l, who mod- 
eled an inhomogeneous ISM including self-gravity, cooling, 
and heating, with a resolution of 40 pc. Low surface density 
gas is stripped easily from any radius, although the overall 
mass loss does not differ much from models with a smooth 
disk. Interpreting the morphology of the remaining gas disk 
had to be done carefully, as runs with different resolutions 
lead to different results. 



7 Magnetic fields 

7.1 Magnetic fields in spiral galaxies 

Evidence for MFs comes mainly from (polarised) radio con- 
tinuum emission (i.e. syncrotron emission of cosmic ray 
electrons spiralling around MFs) and Faraday rotation mea- 
surements. The total intensity of the radio continuum re- 
veals the total strength of the MF in the plane of the sky. 
The polarised intensity and the polarisation angle reveal the 
strength and structure of regular fields in the plane of the 
sky. The observation of MFs is still very difficult, but huge 
improvements in the quality of the data is expected from 
the next generation of radio telescopes, e.g. LOFAR (LOw 
Frequency ARray, and SKA (Square Kilometer Array). 
Typically, spiral galaxies have a regular and a tangled 



MF component (Beck 2005c ) Beck 2007 and references 
therein). The typical total field strength is around 10 /xG. 
The strength of resolved regular fields is a few fj,G. The 
regular field follows the galactic plane and the spiral arms, 
although it is slightly offset. There are also examples of reg- 
ular fields that follow the interarm regions. Occasionally, 
MFs sticking out of the galactic plane are found, usually 
they are associated with supernovae outbursts and galac- 
tic winds. The MF in (gas) spiral arms is mostly tangled. 
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It seems to be connected to the turbulence induced by SF. 
The origins of galactic MFs are not fully understood. A fre- 
quently discussed explanation comes from dynamo models, 



where a seed field is amplified (see e.g. review by Kulsrud 



2005 and references therein). They can explain a toroidal- 



this method can explain succesfully the enhanced polarised 
radio contours found at the upstream edge of RPS galaxies 
HWezgowiec eraL][2007l [Vollmer et aL]|2007] l. This proce- 
dure is applied to a number of RPS candidates (e.g. NGC4522: 
|Vollmeretal.|2006) NGC4654:,Soida et al.|2006ljti| NGC4254: 



Uke regular field to some degree. 



7.2 Magnetic fields in clusters 



In the ICM, evidence for MFs comes from Faraday rota- 
tion measurements, inverse Compton scattering of cosmic 
microwave background photons, and radio halos. Typical 
MF strengths are a few /xG ( [EnBlin et a"L]|2005[ [Feretti & 



|Johnston-Hollitt|2004| and references therein). There is some 
evidence that fields are intermittent and turbulent with auto- 
correlation lengths of a few kpc in cluster centres and a few 
10 kpc in cluster outskirts. Also here, the MFs are probably 
closely linked to ICM turbulence. 

Besides direct measurements, there is plenty of indirect 
evidence for the existence of MFs in the ICM. E.g. bubbles 
rising buoyantly in the ICM seem to be more stable than 
predicted by purely hydrodynamical simulations. Viscosity 
or MFs with coherence lengths larger than the bubble size 
could stabilise them by preventing hydrodynamical instabil- 
ities ( [Reynolds et al.|2005|[Ruszkowski et al.|2007| ). 



7.3 Magnetic fields and ram pressure stripping 

In interacting galaxies, MFs trace regions of gas compres- 
sion, strong shear motions, and enhanced turbulence (Beck, 



2005a I. Such motions can be caused by tidal interactions as 



well as interaction with the ICM. 

Typical features of RPS candidates are asymmetrical dis- 
tributions of the polarised radio intensity, with polarisation 



maxima at the upwind side of the galaxy ( |Wezgowiec et al. 



|2007t I Vollmer et al.|[200T and references therein). Thus, 
the MFs contain information about the velocity components 
in the plane of the sky, which is not accessible otherwise. 
Moreover, distorted MF configurations have been found in 
galaxies that show no other traces of interactions. Thus, 
MFs could be either sensitive even to weak RPs that leave no 
other observable traces, or they may remember distortions 
that happened several 100 Myr ago ( Otmianowska-Mazur| 



Chyzy et al.|2006l|Chyzy|2008 



& Vollmer|2003|fWezgowiec et al.|2007 



The question of MFs and RPS again has two sides: The 
influence of RPS on MFs in spirals is addressed in a series 



of papers based on the method described in Otmianowska 



Mazur & Vollmer (2003 1. First, RPS of a disk galaxy is 



modelled by the sticky particle method of [Vollmer et al.| 
(2001 1. In a manner of post-processing, the initial galaxy 
is given a toriodal MF, which is then evolved by solving the 
induction equation with an MHD grid code using the ve- 
locity fields produced by the sticky particle code. Thus, the 
MF is advected along with the gas, but no effect of the MF 
on the dynamics of the gas is taken into account. Despite 
problems in mapping between the particle and grid code. 



Simulations of the motion of a gas cloud through a mag- 
netised ICM show the effect of draping (Lyutikov 2006; 
[Jones et al.|T996l|Dursi & Pfrommer,2008, Asai et al.,2004i 
[2005|2007||Ruszkowski et al.[2007| l: flie MF is swept up and 
wrapped around the cloud, suggesting that MFs will also 
have an influence on the stripping of disk galaxies. 

A self-consistent approach is still missing. 



8 Intracluster medium properties 

The ICM is an ionised medium, where protons and elec- 
trons have differing mean free paths. Strictly speaking, they 
should be treated as two different fluids. Portnoy et al. ( 1993| l 
investigated the differences between two-fluid and single 
fluid treatments for RPS of elliptical galaxies and found dif- 
ferences in the temperature structure of the wake. 



Due to the strong temperature dependence of the viscos- 
ity and the high ICM temperature, in the unmagnetised case, 
the ICM becomes highly viscous, so that flows on scales of 
galaxies should be viscous. Also based on observations of 
the Perseus cluster, it has been suggested by [Fabian et al. 



( 2003c r), that viscosity may play an important role in dis- 
sipating energy injected by the central AGN. Circumstan- 
tial evidence for the presence of significant ICM viscosity 
is also provided by an examination of the morphology of 
Ha filaments in the Perseus cluster. Several of the filaments 
appear to trace well-defined arcs which argues against the 
presence of strong turbulence in the ICM core, possibly re- 
sulting from the action of viscosity. 

On the other hand, even weak MFs lead to a tiny proton 
gyroradius, which results in a very efficient suppression of 
transport processes like viscosity and thermal conduction. 
If the fields are tangled (Clarke 2004; En BUn et al.|2005 i 



small enough scales, also a fluid description is justified, and 
the macroscopic viscosity may be greatly reduced. 

In summary - the properties of the ICM are poorly con- 
strained ( McNamara & Nulsen|2007 1. In order to study the 
possible impact of a viscous ICM on RPS, Roediger & Briig- 
[gen| ( |2008a) l simulated stripping of a disk galaxy in an ICM 
of various viscosities. Including viscosity requires a very 
small timestep, hence these simulations were done in 2D. 
The stripping radius and amount of lost gas was not affected 
by the viscosity. However, as expected, the more viscous the 
ICM, the smoother the tail of stripped gas. The 2D restric- 
tion prohibits, however, a detailed interpretation of the tail 
structure. 
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9 Enrichment of the intracluster medium by 
ram pressure stripping 

The gas lost from a RP stripped galaxy becomes part of the 
ICM, enriching the ICM with metals. One or more of such 
enrichment processes are expected to take place, because 
the ICM is not primordial but has metalicities of 1/3 to 1/2 
solar Besides RPS, also galactic winds, outflows from ac- 
tive galactic nuclei, and intracluster stars contribute to the 
enrichment. Simulations (see recent reviews by |Schindler &| 
|Diaferio|[2008l [Borgani et aT][2008| find that these enrich- 
ment processes lead to an inhomogeneous metal distribu- 
tion. E.g. RPS preferably pollutes cluster centres, whereas 
galactic outflows can be suppressed by the high ICM pres- 
sure in centres and are more effective in the outskirts. On 
top of the radial variation of the metalUcity in clusters, both, 
observations and simulations (see review by IWerner et al. 



2008 1, show a rich substructure of clumps and filaments. 



Apparently, the lost gas is not mixed with the ICM immedi- 
ately. According to the simulations, no proposed enrichment 
process can produce all observed metals, a combination of 
processes is needed. The upcoming generation of X-ray ob- 
servatories will be able to produce highly resolved metallic- 
ity maps, which will help to constrain and understand the 
various gas removal processes. 



10 Discussion and summary 

Early work as well as more recent simulations agree that in 
many situation the gas removal from disk galaxies by RP 
can be approximated well by the simple analytical G&G 
criterion (Sec. 



2.1 



The estimate works best for galaxies 
moving face-on, and still gives reasonable results for incli- 
nations deviating up to 60° from the face-on case. It fails for 
near-edge on symmetries. Even a time-dependent version of 
the G&G criterion can be applied to galaxies experiencing 
a varying RP along their orbits through a cluster, unless the 
cluster is very compact (Sec. [2]). Hence, RPS seems to be 
able to explain the truncated gas disks found in many clus- 
ter spirals (Sec. [4]l. When RPS removes the gas from the 
outer part of a galactic disk, also SF ceases in the stripped 
parts. This is a plausible idea, and finds general observa- 
tional support. 

Beyond these basics, there is a number of recently stud- 
ied but still unsolved aspects: 

- Several studies (observational and theoretical) suggest 
that RP causes an enhancement in SF activity in the re- 
maining gas disk (Sects. |4] [6|. However, contradicting 
observations exist, and theoretical work needs confirma- 
tion due to the complex nature of SF physics. 

- RPS should produce gas tails behind the stripped galax- 
ies. Some examples are observed, but current simula- 
tions cannot explain the properties of the tail satisfacto- 
rily (Sec. |5]l. 



- The difficulties with the RP tails could arise due to miss- 
ing physics in the models. First extensions towards in- 
cluding MFs (Sec. |7]| and transport processes (Sec. [8| 
are made, but need improvement. Further issues include 
the temperature and ionisation structure inside the tails. 

Advances in these open questions will also assist other 
astrophysical topics: The understanding of the tails of RP 
stripped galaxies will require the inclusion of additional phy- 
sics, most likely ICM transport processes and MFs. In turn, 
success in this question will also help to understand and 
constrain ICM properties, which are relevant for a variety 
of other phenomena, e.g. the interaction between AGN jets 
and the ICM. The question of enhanced SF activity due to 
RP may be used as a test case to check our understanding of 
SF in a special situation. 
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